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INTRODUCTION 

For a number of years, our group has been working on 
generation, manipulation and detection of non-classical 
quantum states of light using techniques and concepts 
from both the continuous variable (cv) and the discrete 
variable (dv) subfields of quantum information science. 
By doing so, we have taken part in a trend - both theoret- 
ical and experimental - of moving away from seeing the 
cv and dv paradigms as separate, competing frameworks 
for quantum information processing and rather look for 
ways of taking advantage of and combining the best prop- 
erties of both. 

In this paper, we review our two most recent results 
in this direction: Entanglement distillation of Gaussian 
states pQ and generation of arbitrary qubit-like squeezed 
states [2 j. These two experiments demonstrate some of 
the possibilities that arise when combining squeezed vac- 
uum or other non-classical cv states with the technique 
of photon subtraction. Before getting to that, we give a 
short (and incomplete) overview of a number of results in 
quantum optics and quantum information that brought 
us to the current state of affairs. 



A BRIEF HISTORY OF BRIDGES 
Motivation 

When optical quantum information processing was de- 
veloped through the 90's, a reasonably strong distinction 
was made between the two different "regimes" of dis- 
crete and continuous variable systems. This distinction 
persists to some extent even today. In fact there is no 
fundamental difference between the two; rather, they are 
complementary ways of describing the states of quantized 
electromagnetic fields [3]. 

Most quantum optics experiments and quantum infor- 
mation schemes consider information contained in dis- 



crete degrees of freedom of an optical field: photon num- 
bers, orthogonal polarization or spatial modes, time bins 
and so on. The research of a number of groups, how- 
ever, looks into states of light that are more conveniently 
expressed in terms of continuous degrees of freedom, typ- 
ically as non-commuting observables like the phase and 
amplitude quadratures or Stokes polarization vectors. A 
typical state representation in dv is the density matrix 
over the available orthogonal modes of the system, while 
cv states are often represented by their Wigner func- 
tions - quasi-probability distributions in the phase space 
of conjugate variables. However, a corresponding den- 
sity matrix can be straightforwardly calculated from a 
Wigner function and vice-versa, so the choice of repre- 
sentation is just one of convenience, clarity and perhaps 
convention. 

The interest in using continuous variables for quan- 
tum information processing originates largely from the 
efficiency of various processes: Squeezed light and coher- 
ent states can be generated deterministically, homodyne 
detection can achieve almost 100% efficiency, and many 
standard operations are deterministic and easy to imple- 
ment in the lab. Furthermore, all of these processes can 
be described in the simple and well understood Gaus- 
sian framework where states are represented solely by 
their means and covariance matrices. Cv quantum in- 
formation processing is not limited to optical systems, 
but have also been explored in atomic systems such as 
spin-polarized ensembles of neutral atoms, and there are 
several promising approaches for efficiently implement- 
ing light-atom interactions [4]. Finally, it should be men- 
tioned that modern optical communications uses quadra- 
ture encoding (cv) and therefore cv approaches to quan- 
tum information may yield better compatibility with ex- 
isting infrastructure. 

Experimentally, one of the driving forces behind opti- 
cal implementations of quantum information primitives 
and protocols has been the process of parametric down- 
conversion. Although there exist protocols for quantum 
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key distribution using only classical coherent states, for 
basically all other purposes a source of non-classical light 
is required, and for that, parametric down- conversion is 
one of the simplest and most efficient methods around. 
By a second-order nonlinear interaction in a crystal or 
waveguide, a single photon from a pump beam is con- 
verted into two daughter photons correlated in a num- 
ber of different degrees of freedom. For weak pump- 
ing, this results in generation of photon pairs, often 
with orthogonal polarization (type II down- conversion) 
such that they are easily separable. For stronger pump- 
ing and degenerate down- converted photons, the single 
output beam has quantum noise fluctuation below the 
standard quantum limit and is said to be in a squeezed 
state. These two different aspects of the down-conversion 
output lend themselves to detection either by photon 
counting/single photon detection or by homodyne de- 
tection, respectively. Photon detection gives a naturally 
discrete measurement outcome, while homodyne detec- 
tion provides a continuous- valued measurement of a given 
quadrature variable. While the photon generation and 
measurement process is probabilistic, squeezing genera- 
tion and homodyne measurement is fully deterministic. 
This was reflected in the first experimental demonstra- 
tions of quantum teleportation by Bouwmeester et al. [5] 
(dv) and Furusawa et al. [6] (cv) - these almost simulta- 
neous papers are perhaps the most prominent examples 
of the distinction between the two paradigms. 

Experiments 

In 2001, Lvovsky et al. [7] presented probably the first 
experiment to truly bridge the dv and cv worlds. They 
generated a single photon state by heralding on detection 
of one of the photons of a parametrically down-converted 
twin pair, and then measured its Wigner function by ho- 
modyne tomography (full state characterization by ho- 
modyne detection of a large ensemble of identically pre- 
pared states) . They essentially measured the continuous- 
variable representation of an intrinsically discrete quan- 
tum state, the one-photon Fock state. Its Wigner func- 
tion showed a prominent dip below zero at the phase 
space origin - a feature that is common to all non- 
Gaussian pure states, but that had not been observed 
in optics until then, since homodyne measurements had 
only been applied to Gaussian states such as squeezed 
vacuum or coherent states. This accomplishment was fol- 
lowed up by Babichev et al. [8 who performed two-mode 
tomography of a dual-rail single photon, and by Zavatta 
et al. [9 who showed the Wigner functions of photon- 
added coherent states, obtained by seeding a coherent 
beam into one output of a non-degenerate parametric 
down-converter. 

A later breakthrough in the merger of dv and 
cv technologies came when Ourjoumtsev et al. [ID] . 



Neergaard-Nielsen et al. [11 J and later Wakui et al. 
[T2] combined squeezed vacuum generation with pho- 
ton detection-heralded homodyne tomography to demon- 
strate small coherent state superposition states, also 
called Schrodinger kittens. Such states of the form 
|^) = | a) — \—ca) will be useful for quantum informa- 
tion processing [T3HTB] but are unfortunately extremely 
difficult to prepare in a deterministic manner. Dakna 
et al. [17] suggested already in 1997 a way to probabilis- 
tically make high-fidelity approximations to these states 
for small coherent amplitudes a, simply by subtracting 
one or more photons from squeezed vacuum. Wenger 
et al. [18 showed the first signs of non-Gaussian statis- 
tics from a state prepared in this way, and this was later 
followed up by the full negative Wigner function state 
preparations [TPUT2] . The photon subtraction operation 
itself has turned out to be a generally useful tool for quan- 
tum information processing and for fundamental tests of 
quantum mechanics, as explored in a series of theoretical 
studies [19H23], 

These experiments originate in the cv world with 
squeezing and homodyne tomography, but incorporate 
photon detection and conditional state preparation from 
the dv toolbox. In the measured output states, their dual 
nature is reflected in the Wigner function which is a mix 
between a single photon state and a squeezed vacuum (it 
is in fact a squeezed photon) . We can therefore no longer 
claim convincingly that the generated states or the exper- 
iments as such belong to the dv or the cv regimes - rather 
they are bridging these two worlds. This trend of dv/cv 
bridging experiments that slowly started with the results 
from Lvovsky and Bellini's labs in the first half of the 
decade has taken off dramatically in the last few years, 
especially since the first kitten generation demonstrations 
in 2006-7. To our knowledge, kitten states have now been 
observed in Paris, Copenhagen, three labs in Tokyo, and 
at NIST. Most recently, photon- number resolving transi- 
tion edge sensors were used to prepare three-photon sub- 
tracted squeezed vacuum [24 and telecom wavelength 
two-photon subtracted squeezed vacuum [25] . 

Other interesting experiments that are worth mention- 
ing are the two-photon Fock state tomography [26 , the 
time-delayed two-photon subtraction for enlarging the 
even kitten state \a) + \—ca) [27], a squeezed kitten state 
prepared by conditional homodyning of a 2-photon Fock 
state [28], remote preparation of entangled kitten states 
[29] , a direct test of the commutation relations, confirmed 
by homodyne tomography [30 , a tunable POVM detec- 
tor combining homodyne detection with photon num- 
ber resolution [31], tomography of photons from a pulse 
pumped OPO [32], tomography of 0-, 1-, 2-photon Fock 
state superpositions [33] and of intermediate single pho- 
ton/squeezed vacuum states [34], noiseless amplification 
of weak coherent states by photon detection [35-38 , pho- 
ton counting-assisted optimal coherent state discrimina- 
tion [39 -41 , and teleportation of a kitten state [42 . Kim 
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FIG. 1. Conceptual setup for generation of single-photon sub- 
tracted squeezed vacuum. 1) Squeezed vacuum preparation. 
2) Photon subtraction. 3) Homodyne tomography. 

[43] has also compiled a nice review of the theoretical 
foundations and experimental progress in this research 
area, while Lvovsky and Raymer [44] gives a comprehen- 
sive review focused on the methods and applications of 
cv state tomography. A more general overview of var- 
ious dv/cv hybrid approaches to quantum information 
processing is provided by van Loock [45] . 

EXPERIMENT BASICS 

Both of the two experiments detailed later are based on 
our setup for generation of photon-subtracted squeezed 
vacuum (kitten state). We will here briefly go through 
some technical details of this experiment, so that we can 
better focus on the conceptual ideas later. 

There are essentially three stages of the kitten gen- 
eration: Resource state preparation (squeezed vacuum), 
state manipulation (photon subtraction) and state char- 
acterization (homodyne tomography). These are illus- 
trated in Figure [I] and described below in turn. For more 
details, refer to [T21I27]. 

Squeezed vacuum 

Our source of non-classical light is a continuously 
pumped optical parametric oscillator (OPO) consisting of 
a 10 mm long periodically poled potassium titanyl phos- 
phate (PPKTP) crystal in a bow-tie type single-ended 
optical cavity. It is pumped by a frequency-doubled 
Ti:Sapph laser with pump wavelength 430 nm. The cav- 



ity has a 580 MHz free spectral range and the FWHM 
bandwidth is 9 MHz, and because of the low internal 
losses, the outcoupling efficiency is as high as 96%, which 
leads to almost pure squeezing at the output. This OPO 
configuration has proven to be able to deliver up to 9 
dB of squeezing [46], but high squeezing levels result in 
a higher degree of mixedness in the observed temporal 
modes because the spectral character of the state be- 
comes "skewed" ; the bandwidth of the squeezed quadra- 
ture is 7 + e, while for the anti-squeezed quadrature it is 
7 — e, where 7 is the OPO bandwidth and e the pump 
parameter |47j [48] . For this reason, we usually generate 
squeezing in the range 2-5 dB. 

In order to control the orientation of the squeezing in 
phase space, we inject a coherent probe beam through 
one of the high-reflectivity mirrors of the OPO. This 
beam then co-propagates with the squeezed vacuum and 
can be used for phase locking as well as locking of the 
filter cavities in the photon subtraction channels. For 
phase locking the squeezing orientation, we monitor the 
transmitted beam on a standard DC detector. The beam 
has undergone classical phase-dependent parametric am- 
plification in the OPO, so by fixing the amount of am- 
plification, the probe beam will be locked in phase with 
the squeezed vacuum. By later locking the probe beam 
to the local oscillator of the homodyne detector (which 
determines the phase space orientation), we can fix the 
squeezing orientation. 

The probe beam is much stronger than the squeezed 
vacuum, and since the photon detector (APD) is fre- 
quency undiscriminating and cannot distinguish between 
the narrow-band probe and the wide-band squeezing, we 
need to avoid the probe beam when measuring the kit- 
tens. Therefore we introduced a chopping cycle of the 
probe beam where it is on for 20 [is and off for the next 80 
/is, chopped by two double-pass AOMs. While the probe 
is on we do the phase and cavity locking and shield the 
APD. While it is off, we uncover the APD and measure 
the conditional homodyne signals. In the meantime, the 
locks are kept passively stable. 

Photon subtraction 

The squeezed light impinges on a 5% reflection beam- 
splitter, where the main transmitted part goes directly 
to the homodyne detector for output state analysis. The 
small reflected part is directed to an avalanche photo 
diode (APD), whose click signal is connected to a digital 
oscilloscope to act as a conditioning trigger of the homo- 
dyne data acquisition. The squeezed light does not con- 
tain many photons, only 5% is tapped for trigger detec- 
tion, and the total transmission and detection efficiency 
of the trigger channel is only about 10%, so the pho- 
ton count rate is rather low - around 10,000 counts/sec 
for 3 dB squeezing. But whenever a photon is detected, 



we know that it is missing from the correlated part of 
the main signal beam, which is therefore projected to 
a photon-subtracted squeezed vacuum state. The low 
efficiency of the trigger channel only influences the suc- 
cess rate of the experiment; the impact on the generated 
states is almost negligible. 

An APD and a homodyne detector have completely 
different spectral and temporal mode resolutions, so the 
photon detected by the APD is not necessarily correlated 
with the mode seen by the homodyne detector. That 
is true for both a pulsed laser setup [10 and for a cw 
setup [TTJ [12] , but the problem is manifested in differ- 
ent ways and requires different solutions. In the cw case, 
the down- conversion process happens over a nanometer- 
broad spectral range, but only at those frequencies res- 
onant with the OPO cavity. The APD can detect all 
the down-converted photons, but the homodyne detec- 
tor is limited to a bandwidth of around 50 MHz around 
the carrier frequency of the local oscillator, derived from 
the main 860 nm laser. Therefore we have to eliminate 
all the non-degenerate OPO spectral cavity modes from 
the light reaching the APD. We do that by passing the 
light through two subsequent Fabry-Perot cavities, both 
resonant with the degenerate OPO mode but with differ- 
ent, rather wide free spectral ranges. This arrangement 
works very well and seemingly filters away all unwanted 
photons. 

Homodyne tomography 

To measure the Wigner functions of the generated non- 
classical light states, we interfere them with a strong lo- 
cal oscillator (LO) on a 50:50 beam splitter, detect the 
two output beams on fast PIN photo diodes and mon- 
itor the amplified difference photocurrent on the oscil- 
loscope triggered by the APD signal. For every trigger 
click, we record a trace of about 1 /is length around the 
trigger time at a sample rate faster than twice the detec- 
tor bandwidth. Later, in post-processing, we perform a 
temporal mode selection by weighting each trace with a 
mode function 

fit) = Ke-^-^ - 7 e-* :|t - £o1 , (1) 

where 7 is the OPO HWHM bandwidth, k the filter 
cavity bandwidth, and to the time of photon detection. 
This mode function corresponds to the temporal extent 
of the correlations between the detected photon and the 
photon-subtracted main beam, so it is close to the opti- 
mal temporal filter for extracting the conditionally pre- 
pared state with high efficiency. 

The homodyne detection process described here results 
in a single outcome of the quadrature measurement at 
the given LO phase. To reconstruct the full state, we 
repeat the measurement many thousand times for iden- 
tically prepared states, with the LO phase locked to one 
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FIG. 2. Example of homodyne tomography of a single-photon 
subtracted squeezed vacuum state. On the left is a series 
of homodyne measurements taken at 12 different fixed phase 
space angles. Notice the relative lack of observations around 
the origin - this is reflected in the negative dip of the re- 
constructed Wigner function in the right panel. The Wigner 
function has been rotated in phase space to have squeezing 
along the p-axis. 




FIG. 3. The Gaussian toolbox. A basic continuous- variable 
experiment involving the essential Gaussian operations. For 
a number of important quantum information tasks, this set 
of Gaussian tools is not sufficient and must be supplemented 
by a non-Gaussian tool such as a single photon detector. 

of 6 or 12 evenly distributed values between and tt. 
Finally, we process the total set of quadrature measure- 
ments with a maximum likelihood algorithm [44j EH] to 
estimate the density matrix of the generated state. An 
example of a full tomographic measurement along with 
the reconstructed state (Wigner function calculated from 
the density matrix) is presented in Figure |2j 

DISTILLATION OF GAUSSIAN 
ENTANGLEMENT 

One of the great advantages of cv quantum informa- 
tion processing is that it is essentially constructed from 
deterministic, Gaussian operations, that is, operations 
with Hamiltonians that are at most quadratic in the field 
operators. All of the most commonly used techniques in 



5 



the lab such as squeezing, phase shifting, beam split- 
ting, displacement and homodyne detection (Figure [3| 
are Gaussian operations, which means that an input state 
with Gaussian Wigner function (or characteristic func- 
tion) remains Gaussian after the transformation. The 
theory for Gaussian states, entanglement, operations and 
so on is very well developed and provides a powerful 
framework for various quantum information processing 
tasks [50j ISI] • Unfortunately, there are some significant 
exceptions to what is possible with Gaussian resources. 
Bartlett et al. [52] showed that a quantum computer us- 
ing only Gaussian resources can be efficiently simulated 
by a classical computer. Therefore non- Gaussian opera- 
tions are needed to achieve a quantum-mechanical speed- 
up of computing. Another problem is long-distance dis- 
tribution of entanglement which is needed for quantum 
communication protocols. In order to combat propaga- 
tion losses, the entanglement must be reinforced by quan- 
tum repeater stations [53], incorporating entanglement 
swapping and distillation [53] , However, it was shown 
that a Gaussian entangled state cannot be distilled by 
Gaussian local operations and classical communication 
[55H57] . Since Gaussian entanglement is the predomi- 
nant type of cv entanglement (such as two-mode squeezed 
EPR states), non-Gaussian operations are required for 
distribution over long distances. 

Hage et al. [58, 59 and Dong et al. [60 performed en- 
tanglement distillation using only conditional homodyne 
detection, but in those cases, the Gaussian entangled 
beams had been exposed to non-Gaussian noise. Then 
it is possible to recover the original amount of Gaussian 
entanglement by Gaussian operations such as homodyne 
detection. However, the probably most common and crit- 
ical transmission noise comes from linear loss, which is a 
Gaussian transformation. Ourjoumtsev et al. [61] demon- 
strated that a non-local (collective) photon subtraction 
can increase entanglement that is still Gaussian. Based 
on this demonstration and on the realization by Opatrny 
et al. [19] that photon subtraction can increase entangle- 
ment and lead to better teleportation, we performed the 
first entanglement distillation Q of a Gaussian two- mode 
state pp. 

Experiment 

The basic concept of the distillation protocol, illus- 
trated in Figure [4j is to prepare a series of entangled 



We note that there is some confusion about terminology in the 
literature about entanglement distillation: By some accounts, 
our result should rather be classified as entanglement concentra- 
tion, while the term distillation is in line with other sources - 
notably the conceptually similar procrustean method for qubit 
distillation of Kwiat et al, Nature 409, 1014 (2001). 




FIG. 4. Basic concept of entanglement distillation. Weak 
entanglement is enhanced by probabilistic local operations 
(orange boxes) followed up by classical communication (blue 
wires) to select those entangled pairs that have been success- 
fully distilled. 




FIG. 5. Entanglement distillation experimental setup, in this 
case for simultaneous photon subtractions at both Alice and 
Bob. The Wigner functions illustrate the measured and re- 
constructed states of the A-B mode. The A+B mode is in all 
cases vacuum. 

states shared by Alice and Bob, let them perform local 
operations on their part and finally let them communi- 
cate with each other classically to determine which states 
have successfully had their entanglement increased and 
therefore can be kept for further processing. This pro- 
cess sacrifices the rate of shared states for the amount of 
entanglement. 

In our experiment (Figure [5| we generate a two- mode 
entangled state simply by splitting a squeezed vacuum 
of variable squeezing level on a 50:50 beamsplitter. This 
creates a state with rather weak entanglement, equivalent 
to that of a standard two-mode squeezed (EPR) state 
with half as much squeezing. For Alice and Bob's local 
operations, we use the photon subtraction technique de- 
scribed in the previous section, with a separate channel of 
frequency filters and APD for both. A click of the APD 
detectors signals a successful distillation event. If Alice 
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FIG. 6. Full two- mode homodyne tomography of a single- 
party distilled state. On the left are the joint quadrature 
distributions for all 36 combinations of Alice's and Bob's LO 
phases. From these data the full density matrix (right) was 
reconstructed. When combining the density matrix elements 
in groups with the same total photon number, the overall 
block structure resembles a single mode photon-subtracted 
squeezed vacuum. 



and Bob were to proceed with an application of their 
shared entanglement, they would have to classically tell 
each other at what times they got a detector click. We do 
not go that far here - instead we simply characterize the 
resulting quantum state. To quantify the entanglement 
before and after the distillation, we reconstruct the joint 
two-mode density matrix by homodyne detection on both 
output modes for both the undistilled state and for the 
state with a photon subtracted simultaneously from each 
of the two modes. It turns out that even subtracting a 
photon from just Alice or Bob also increases their shared 
entanglement, so we did a third series of measurements 
of this scenario. The photon subtraction beamsplitter 
reflectivities were set at 0% for the initial state measure- 
ment, at 5% for the single photon detection series and 
at 10% for the two simultaneous photon detection series 
(as a compromise between success rate and entanglement 
gain). 

It is possible to do a complete two-mode tomographic 
measurement by independently varying the LO phases of 
Alice's and Bob's homodyne detectors. Figure [6] shows 
an example of a density matrix reconstructed from such 
a measurement, where 36 LO phase combinations (0-0, 
O-tt/6, 0-5tt/6, tt/6-0, 7r/6-57r/6, 57r/6-57r/6) 
were used. However, this process is very time-consuming, 
and due to our particular kind of initial entangled state, 
the output state can be separated into a vacuum state 
in the xa + xb,Pa + Pb variables and a state in the 
%a — %b,Pa — Pb variables that corresponds to a photon- 
subtracted squeezed vacuum. This fact allows us to do 
the two-mode tomography using only the 6 phase com- 
binations where the two LO's are in phase, a strategy 
that is experimentally justified by confirming that the 
virtual mode is indeed in a vacuum state. For de- 
tails, refer to the Supplementary Information of [I]. In 
broad terms, the argument relies on the fact that the 



experimental setup is fully equivalent to a setup where 
the photon subtraction (by 1 or 2 APDs) is done before 
the 50:50 beam splitting. That same fact shows that for 
the single-party distillation operation, the joint output 
state is independent on whether Alice or Bob detected 
the photon. 



Results 

The result of the distillation is shown in Figure [7^i. 
The entanglement measure we use is the logarithmic neg- 
ativity [62] which is easily computable from the density 
matrix. The graph shows that for a wide range of initial 
squeezing levels, the entanglement gain from the distilla- 
tion protocol is significant for both single-party and dual- 
party operations. Perhaps surprisingly, the single-party 
distillation mostly works better than the dual-party. This 
is a result of the particular kind of entangled initial state 
we use. It can be understood by realizing that for weak 
squeezing levels, the shared state after a single photon 
subtraction is essentially a split single photon which has 
a full e-bit of entanglement, whereas for two photon sub- 
traction (or zero) , the shared state is mostly vacuum with 
a small contribution of photon pairs and therefore with 
entanglement close to zero. However, the effect of the 
single-party distillation operation is not limited to this 
particular scenario - the pure-state calculations in Fig- 
ure [7]3 show that also for the more general two-mode 
squeezed vacuum state there is a marked increase in en- 
tanglement for both single-party and dual-party photon 
subtractions. 

For a complete Gaussian input to Gaussian output 
distillation protocol, the procedure demonstrated here 
should be expanded to include several similarly prepared 
states and then be combined with an iterative Gaussifica- 
tion procedure [20]. However, even after just this single 
step the resulting distilled state is immediately useful for 
e.g. improvement of teleportation fidelity [I9j EH [64] , at 
least for the dual-party photon subtraction. This can be 
seen from Figure [7)3, showing that the two-click distilled 
state not only has higher entanglement, but also a higher 
degree of squeezing than the input. 



CONTINUOUS VARIABLE QUBIT STATES 

The discovery by Knill et al [65] that quantum com- 
putation is possible using only linear optics in combi- 
nation with single photon sources and photon detectors 
gave birth to a wealth of research into linear optics quan- 
tum computation (LOQC). Specifically, a number of pro- 
posals suggested to use coherent states as the basis for 
circuit-style LOQC [66-68 ; that is, to encode logical 
qubits into the physically higher-dimensional coherent 
states \a) and |— a). It was shown that such schemes 
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FIG. 7. Entanglement distillation, a) Logarithmic negativity for the undistilled, single-party (1 click) distilled and dual- 
party distilled states for various input squeezing levels. Dashed curves are calculated from a model taking all experimental 
parameters into account. The 3D plots are the reconstructed Wigner functions W(xa—xb,Pa ~Pb) for 3.2 dB initial squeezing, 
b) pa ~Pb quadrature variance for the undistilled and the dual-party distilled states. Because of losses and other imperfections, 
the undistilled output squeezing is lower than the initial level, whereas the distilled state is actually more strongly squeezed 
than the initial state, c) Theoretically calculated entanglement negativity for pure states of both 50:50 split squeezed vacuum 
as used in the experiment and of standard two- mode squeezed vacuum which initially has twice as much entanglement. Both 
single-party and dual-party photon subtraction are effective for both classes of entangled state. 



could allow efficient computation with only simple linear 
operations in-line, with the hard operations put off-line 
to the preparation of the resource states. For large am- 
plitudes a the two states are basically orthogonal, but 
for small a they have a finite overlap - nonetheless, this 
scheme is capable of fault-tolerant quantum computing 
with certain advantages (as well as disadvantages) over 
traditional single photon-encoded schemes [16] . One par- 
ticular advantage is - in contrast to the usual two qubit 
case - that all four Bell states can be unambiguously dis- 
criminated, leading to higher efficiency for e.g. quantum 
teleportation [23]. It is also noteworthy that for long- 
distance optical communication, it has been shown that 
simple coherent states attain the channel capacity, but 
the receiver would have to employ collective decoding 
with a quantum computer for coherent state signals to 
extract the maximum information [69] . 

To realize a coherent state-based computing scheme 
as described here, we need access to resources of arbi- 
trary qubits a\a) + b\—a). The diagonal states with 
a = 1,6 = ±1 (unnormalized) were already approxi- 
mately realized in the various kitten state experiments 
mentioned in the introduction, and of course the two ba- 
sis states are trivial, but all other qubit states still remain 
to be demonstrated with high fidelities. We suggested a 



way to accomplish the generation of these complex co- 
herent state superpositions using two-photon subtraction 
of squeezed vacuum with an added displacement [70] . 
The displacement before the second photon detector can 
be adjusted in phase and amplitude, according to which 
the conditional output state is prepared in any arbitrary 
superposition of one-photon and two-photon subtracted 
states (basically equivalent to a rotation of the coherent 
state qubit space) as long as a is small. 



In our experimental study [2 , we demonstrated this 
method of displaced photon subtraction, although we 
used only a single photon detector. Our generated states 
are therefore superpositions of squeezed vacuum and 
single-photon subtracted squeezed vacuum (equivalent to 
a squeezed photon). The concept is illustrated in Figure 
[8] In a standard kitten state generation setup we in- 
sert a highly reflective beamsplitter in the trigger beam 
just before the APD. On that beamsplitter we inject a 
strong coherent beam with variable phase and amplitude 
which acts to displace the state of the trigger beam in 
phase space. Now the photon detection on the APD be- 
comes ambiguous; there are two possible, indistinguish- 
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FIG. 8. Superposition state generation by displaced photon 
subtraction, a) Standard odd kitten state factory, b) With 
the trigger beam from the squeezing blocked, all photon clicks 
are from the uncorrelated auxiliary beam. The output state is 
then just squeezed vacuum, c) With trigger beam and aux- 
iliary displacement beam overlapped on a highly reflective 
beamsplitter, the two possible sources of the photon detec- 
tion are indistinguishable and the output is projected to a 
superposition state determined by the relative strength and 
phase of the two beams. 

able sources of the detected photon[^J so the total effect of 
a click on the trigger beam is a projection onto /3|0) + |1) 
and therefore the conditional output state ends up in a 
superposition of single-photon subtracted squeezed vac- 
uum and normal squeezed vacuum. The exact parame- 
ters of the superposition depend on the phase and ampli- 
tude of the displacement beam. If the displacement beam 
is strong (weak), the output will be close to a squeezed 
vacuum (photon). 

To illustrate the potential of this procedure, we carried 
out a large number of state preparation and character- 
izations for different displacement parameters. Because 
we restricted ourselves to single-photon subtraction, we 
could perform detailed homodyne tomography (360,000 
samples at 12 fixed phases) for all these states within a 
reasonable time span. Figure [9] shows the density ma- 



2 The bandwidths of the squeezed vacuum and the coherent dis- 
placement beam are quite different, but the very wide integrated 
bandwidth of the APD erases this source of distinguishability. 




FIG. 9. Fock state basis density matrices of representative 
reconstructed output states. From left to right: Squeezed 
vacuum, squeezed photon, a state generated with a displace- 
ment beam of 1/3 the intensity (at the APD) of the squeezed 
trigger beam and with 180° relative phase, and finally a simi- 
lar state but with 90° relative phase. The upper matrices are 
real parts, lower matrices are imaginary. 




FIG. 10. Squeezed Fock state basis density matrices of the 
same states as in Figure [9] 

trices of just a few of these states: squeezed vacuum, 
squeezed photon, and two states roughly halfway in- 
between where the trigger and displacement beams were 
at 180° and 90° relative phases. No correction for exper- 
imental inefficiencies were performed (and the same goes 
for the previous distillation results). We have established 
a detailed model of the full experiment that fits the mea- 
sured outcomes very well without any free parameters, 
showing that we have a very good understanding of the 
individual components of the setup and providing guides 
to how to improve the results [48J. 

This experiment truly bridges the dv and cv worlds. 
Not only do we use experimental methods from both, 
but the output states seen as a whole can also be consid- 
ered as qubits of squeezed states, being superpositions of 
the two orthogonal states S(r)\0) and S(r)\ 1). To clearly 
see this, we have transformed the density matrices to a 
different basis in Figure [ToJ namely the squeezed Fock 
state basis |a\) = S(r)\k) - since the squeezing opera- 
tor is unitary, this is a proper transformation where the 
new basis is orthonormal and complete. In this squeezed 
basis, the states are essentially confined to the lowest 2- 
dimensional subspace {|</>o)> with a small contribu- 

tion from higher-order states mostly due to experimental 
inefficiencies. From these density matrices, we can calcu- 
late the parameters for each state when interpreted as a 
squeezed qubit 

|V)=cos^(r)|0)+e^sin^(r)|l). (2) 
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For the two non-trivial states of Figures [9| [10] the param- 
eters are = 100° and (j> = 0°,90°, respectively. The 
states are not pure - they range between 0.56 to 0.88 in 
purity - and therefore the squeezed vacuum and squeezed 
photon are also not perfectly orthogonal, with an overlap 
of 0.26. The generated qubit states are therefore not di- 
rectly usable for quantum information purposes as such, 
but they still serve as a clear demonstration of the po- 
tential of the displaced photon subtraction method. In 
fact, it is also applicable to genuine cat states, since the 
photon subtraction operation a turns the even cat state 
| a) + |— a) into the odd \a) — \— a) and vice versa, and 
the displaced photon subtraction can therefore generate 
an arbitrary superposition of these two states. It was 
also pointed out by Marek and Fiurasek [71] that the op- 
eration corresponds to a single-mode phase gate in the 
coherent state basis, and that in more complex settings, 
it can be used to implement the Hadamard gate and a 
two-mode phase gate, which just serves to emphasize its 
generally useful quality. 

To finish this section and to emphasize the dv-cv link, 



we present in Figure [TT] the full set of generated qubit 
states, represented by their reconstructed Wigner func- 
tions (cv) and inserted at their appropriate positions 
in a Bloch/Poincare sphere (dv) according to the qubit 
parameters extracted from their squeezed basis density 
matrices. To see the detailed parameters and close-up 
Wigner functions of each state, please go to our group's 
website where an interactive Java applet is available [72] . 



CONCLUSION 

We have presented two different examples of what can 
be accomplished with the photon subtraction operation 
on non-classical states of light. From a practical appli- 
cations point of view, it enables a number of different 
quantum information protocols in continuous variable or 
hybrid regimes. But it also highlights that the division 
between the discrete variable and continuous variable 
worlds is not as large as it was previously made out to 
be, perhaps to the extent that it hardly makes sense to 
distinguish them. Judged from the amount of activity on 
these kinds of hybrid processes and schemes, both on the 
theoretical and experimental side, we expect to see many 
more breakthroughs in the near future, ultimately lead- 
ing to practically applicable quantum information tech- 
nologies. 
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